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Abstract 


Although  innumerable  different  ionic  liquids  are  possible,  even  basic  physical-property  data, 
such  as  the  density  and  melting  point,  exist  only  for  relatively  few.  Derivation  of  melting  point 
QSPRs  (Quantitative  Structure-Property  Relationships)  for  energetic  ionic  liquids  would 
therefore  greatly  aid  in  the  molecular  design  of  new  compounds.  A  new  class  of  ionic  liquids, 
based  on  l-substituted-4-amino-l,2,4-triazolium  bromide  and  nitrate  salts,  were  recently 
synthesized  and  their  melting  points  and  densities  measured.  We  optimized  the  molecular 
geometries  of  the  cations  of  the  ionic  liquids  using  ab  initio  quantum  chemical  methods. 
Melting-point  QSPRs  were  then  derived  from  molecular  orbital,  thermodynamic,  and 
electrostatic  descriptors.  Good  correlations  with  the  experimental  data  were  found.  The 
correlation  coefficients  for  three-parameter  melting-point  QSPRs  and  for  one -parameter  density 
QSPRs  exceed  0.9.  Although  some  of  the  descriptors  that  appear  in  our  QSPRs  were  designed  to 
describe  chemical  reactions,  we  infer  that  they  serve  in  this  study  to  quantify  interactions 
between  the  cation  and  anion. 
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Introduction 


Hydrazine  and  its  methylated  derivatives  are  powerful  reducing  agents  with  a  wide  range  of 
uses,  including  aerospace  fuels.1  The  introduction  of  hydrazine,  monomethyl  hydrazine,  and  1,1- 
dimethyl  hydrazine  grew  out  of  a  need  for  high-energy,  noncryogenic,  liquid  fuels  that  can  be 
used  alone  or  mixed  with  other  components.  The  toxicity  of  hydrazine  fuels  is  a  substantial 
operational  concern  to  the  U.S.  Department  of  Defense  as  well  as  to  the  aerospace  industry.  The 
two  most  likely  occupational  exposure  routes  are  inhalation  and  skin  exposure."  We  have 
previously  reported  quantitative  structure-activity  relationships  for  toxicity  prediction  of  a  series 
of  high-energy  compounds. 

Ionic  liquids  have  melting  points  that  fall  below  room  temperature  and  have  generated 
considerable  interest  as  ‘green’  alternatives  in  current  industrial  applications  because  of  their 
exceptional  solvating  ability  and  extremely  low  vapor  pressures. 4’5’6,7,8  Energetic  ionic  liquids 
should  be  inherently  less  toxic  than  hydrazine-based  fuels  because  the  exposure  routes  are 

1  o 

greatly  diminished  due  to  their  low  vapor  pressures.  Although  approximately  10  different  ionic 
liquids  are  possible,9  even  basic  physical-property  data,  such  as  the  density  and  melting  point, 
exist  for  relatively  few.  Derivation  of  melting-point  QSPRs  (Quantitative  Structure-Property 
Relationships  -  models  that  relate  chemical  structure  to  physical  properties)  for  energetic  ionic 
liquids  would  therefore  greatly  aid  in  the  molecular  design  of  new  families  of  energetic  ionic 
liquids.  QSPRs  for  the  prediction  of  densities  of  energetic  ionic  liquids,  which  are  correlated  to 
their  energy  densities,  would  also  be  beneficial. 

Melting  points  for  organic  molecules  depend  in  general  on  the  arrangement  of  atoms  in  the 
crystal  lattice  as  well  as  on  pairwise  interactions.  Difficulty  in  predicting  melting  points  arises, 
in  part,  from  the  fact  that  these  two  properties  are  interdependent.  Currently,  the  descriptors 
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available  for  derivation  of  melting-point  QSPRs  were  not  designed  for  the  description  of 
condensed  media.10  However,  additional  factors  that  affect  the  strength  of  a  crystal  lattice  and, 
therefore,  the  melting  point  are  molecular  symmetry,11  the  molecule’s  conformational  degrees  of 

freedom11  (both  of  which  can  be  accounted  for  by  descriptors  obtained  from  gas-phase  quantum 
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mechanical  calculations),  and  the  molecular  motion  in  crystals. 

Inorganic  ionic  compounds  typically  display  high  melting  points  due  to  strong  electrostatic 
intennolecular  forces.  The  intennolecular  forces  affecting  the  melting  point  of  organic 
compounds  are  much  weaker,  e.g.,  hydrogen  bonding  and  dipole-dipole  interactions.  The  most 
accurate  (RMS  error  =  0.5°C)  melting-point  QSPR  for  organic  compounds,  which  was  derived 
for  24  unbranched  alkanes,  is  composed  entirely  of  topological  indices.13  However,  the  RMS 
error  for  the  QSPR  predicting  the  melting  point  of  56  normal  and  branched  alkanes,  which  also 
used  connectivity  indices,  is  much  larger  (23.8°C).14  In  a  study  of  303  nonnal  and  branched 
substituted  alkanes,  intennolecular  forces  provided  descriptors  for  an  11 -term  QSPR  with  an 
RMS  enor  of  16.4°C.15  In  developing  QSPRs  for  melting  points  of  substituted  benzenes, 
Katritzky,  et  ah, 16  found  that  a  hydrogen-bonding  descriptor  together  with  other  quantum 
chemical  descriptors  led  to  correlations  much  improved  over  QSPRs  derived  with  only 
topological,  geometrical,  and  other  traditional  descriptors. 

Descriptors  that  characterize  weak  intennolecular  interactions  are  expected  to  provide  accurate 
melting-point  QSPRs  for  ionic  liquids  because  they  are  weakly  coordinating.  Melting-point 
QSARs  derived  by  Katritzky  and  coworkers  for  ionic  liquids  and  ionic-liquid  analogues  (based 
on  pyridinium  bromides17  and  on  imidazolium  and  benzimidazolium  bromides9)  are  composed  of 

IS. 

information  content  indices  -  descriptors  characterizing  the  cations’  size  and  electrostatics  - 
and  the  average  nucleophilic  reaction  index,  i.e.,  the  propensity  of  the  cations  to  react  with 
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nucleophiles.19  The  resulting  melting-point  QSPRs  have  r2  values  ranging  from  0.690-0.788, 
except  for  the  QSPR  derived  for  a  subset  of  18  substituted  imidazolium  bromides,  which  has  an 
r2  value  of  0.943. 9,17 

20  2 1 

Liquid-density  QSPRs  have  been  previously  derived  for  alkanes  in  two  studies."  ’  In  the  first 
study,  structural  descriptors  computed  from  molecular  graphs  were  used  to  derive  physical 
properties,  including  liquid  density,  for  134  alkanes.  Using  the  same  dataset,  better  physical- 
property  QSPRs  were  derived  using  descriptors  based  on  graphs  of  atomic  orbitals."  In  neither 
study  were  physical  interpretations  of  the  liquid-density  QSPRs  attempted. 

A  new  class  of  energetic  ionic  liquids,  based  on  l-substituted-4-amino-l,2,4-triazolium 
bromides  and  nitrates  (see  Figure  1),  was  recently  synthesized  and  their  melting  points  measured 
(see  Table  1).  "  Some  of  the  data  in  Table  1  actually  correspond  to  glass  transitions  because  the 
glass  forms  much  more  readily  than  the  crystal.  In  these  cases,  we  employ  glass  transition 
temperatures  in  lieu  of  melting  points  because  we  are  interested  in  the  temperature  at  which  the 
salt  ceases  to  be  a  liquid.  As  reported  by  Drake,  et  al.,  "  these  new  ionic  liquids  are  stable  at 
room  temperature  and  show  no  signs  of  decomposition  after  more  than  a  year  at  ambient 
temperatures.  The  nitrate  salts  generally  have  lower  melting  points  than  the  bromide  salts.  For 
the  //-alkyl-substituted  compounds,  the  melting  points  generally  increase  as  the  length  of  the 
alkyl  group  increases  but  this  increase  is  not  monotonic;  compounds  with  an  even  number  of 
substituent  carbons  generally  display  lower  melting  points  that  those  with  an  odd  number."" 
Densities  of  the  bromide  salts  were  also  reported  by  Drake,  et  al.,  "  cf.  Table  1.  We  endeavor  in 
this  paper  to  derive  melting-point  and  density  QSPRs  for  these  new  energetic  ionic  liquids  and 
interpret  them  in  tenns  of  molecular  interactions. 
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Methods 


All  ab  initio  quantum  mechanical  calculations  were  performed  with  Gaussian98,  version  A.9.2  ' 
The  optimal  geometry  for  each  cation  (Table  1)  was  obtained  using  the  restricted  Hartree-Fock 
method  (RHF)  and  the  6-3 1G**  basis  set  (6-3 1G**  is  a  split-valence  basis  set  with  polarization 
functions  on  all  atoms).  Vibrational  frequencies  were  calculated  for  all  cations  using  the  same 
level  of  theory  as  the  optimizations. 

CODESSA  version  2.6 124  was  employed  to  derive  correlations  between  each  descriptor  and 
the  melting-point  data,  to  derive  the  QSPRs,  and  to  calculate  the  statistics  for  the  QSPRs, 
including  the  correlation  coefficient,  r,  the  Fisher  significance  parameter,  F,  and  q,  the  cross- 
validated  correlation  coefficient  calculated  using  a  leave-one-out  method.  CODESSA’s  best 
multilinear  regression  method  and  heuristic  method  were  used  to  obtain  three-parameter  melting- 
point  QSPRs.  The  heuristic  method  was  used  to  obtain  one -parameter  density  QSPRs.  Both 
methods  limit  the  collinearity  of  the  selected  descriptors  and  utilize  statistical  significance  to 
derive  QSPRs. 

CODESSA  categorizes  descriptors  as  constitutional,  topological,  geometrical,  electrostatic, 
quantum  mechanical,  and  thermodynamic,  although  both  electrostatic  and  thennodynamic 
descriptors  are  obtained  from  quantum  mechanics  calculations.  Constitutional,  topological, 
geometrical  descriptors  were  seldom  found  to  be  among  the  best  descriptors  and  rarely  improved 
QSPRs  that  were  derived  without  them.  In  this  work,  we  therefore  employed  only  electrostatic, 
quantum  mechanical  and  thermodynamic  descriptors  to  derive  QSPRs.  The  descriptor  pool 
consisted  of  about  250  descriptors.  Unless  specified  otherwise,  the  partial  atomic  charges  that 
were  employed  were  calculated  using  the  Natural  Bond  Orbital  method,  as  implemented  in 
Gaussian’ 03. 23 
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Results  and  Discussion 


Even  though  they  pertain  to  the  gas-phase,  optimal  molecular  geometries  of  the  cations 
compare  well  to  the  crystal  structures  of  the  ethyl,  //-propyl,  //-hexyl,  and  /z-heptyl  bromide  salts, 
which  were  the  only  bromide  salt  crystal  structures  published,-'  and  to  the  isopropyl  and 
methylcyclopropyl  nitrate  salts,  which  were  the  only  two  nitrate  salt  crystal  structrures 
published.-  As  seen  in  Table  2,  calculated  bonds  lengths  are  in  good  agreement  to  those  found 
in  the  crystal  structures.  Although  only  bonds  common  to  all  cations  are  included  in  Table  2, 
experimental  and  theoretical  bond  lengths  in  the  substituents  are  also  in  good  agreement.  A 
linear  fit  of  the  bond  lengths  in  Table  2  has  an  r  value  of  0.964.  The  average  RMS  deviation  for 
the  bond  lengths  in  Table  2  is  0.0163  A. 

Good  correlations  with  the  experimental  melting-poing  data  were  found.  The  following  QSPR 
was  derived  for  the  melting  points  (Tm)  of  the  bromide  salts: 

Tm  =  -262.  -  6.91  x  105  NRInh2  +  47.4  HACA2  -  136./ELUMO 

(1) 

N=13,  r=0.914,  F=31.9,  s2=26.6,  q2=0.784 

where  NRfym  is  the  nucleophilic  reactivity  index19  for  the  amine  nitrogen  (see  Figure  1), 
HACA2  is  the  area-weighted  surface  charge  of  hydrogen-bond  acceptor  atoms,  and  ELumo  is 
the  energy  of  the  lowest  unoccupied  molecular  orbital.  Values  for  all  melting-point  descriptors 
are  presented  in  Table  3.  Correlations  between  descriptors  appearing  in  Eqn.  1  are  poor,  as  seen 
in  Table  4a.  Although  NRI,  which  is  essentially  a  Fukui  function,  is  implicitly  designed  to 
describe  chemical  reactions,19  we  infer  that  it  serves  here  to  quantify  interactions  between  the 
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cation  and  anion.  Similarly,  reciprocal  ELUmo,  as  we  have  shown  previously,28  can  be  interpreted 
within  the  context  of  frontier  molecular  orbital  theory-  to  describe  the  stabilization  due  to 
mixing  of  occupied  MOs  of  a  molecule  with  the  LUMO  of  the  cation.  Therefore,  reciprocal 
Elumo  may  also  serve  here  to  quantify  interactions  between  the  cation  and  anion.  It  is  also 
plausible  that  the  descriptors  in  Equation  1  account  for  interactions  between  cations.  The  RMS 
error  in  Equation  1  is  about  5°C. 

The  following  QSPR  was  derived  for  the  Tm  of  the  nitrate  salts: 


Tm  =  -284.  -  214.  HDCAlz  -3.94  x  104  NRImm,c  +  3.16  x  103  FHDCA 


(2) 


N=13,  r2=0.933,  F=41.5,  s2=198,  q2=0.872 


26 

Both  HDCAlz  and  FHDCA  are  measures  of  the  hydrogen-bond  donating  ability  of  the  cation 
and  NRImin  r  is  the  minimum  NRI  for  a  carbon  atom.  HDCA1Z  is  the  sum  of  the  hydrogen-bond 
donor  surface  area  weighted  by  each  donor’s  partial  charge  and  FHDCA  is  HDCA1  divided  by 
the  total  molecular  surface  area,  i.e.,  the  fractional  HDCA1.  Both  HDCAlz  and  FHDCA1  are 
interpreted  as  hydrogen-bond  donating  ability.  Equation  2  is  similar  to  Equation  1  in  that  two 
tenns  include  a  NRI  and  a  hydrogen-bond  descriptor.  As  in  Equation  1 ,  we  infer  that  NRImm  c 
and  the  hydrogen-bond  descriptors  in  Equation  2  serve  here  to  quantify  interactions  between  the 
cation  and  anion  or  between  cations.  See  Table  3  for  values  for  all  melting-point  descriptors. 
Correlations  between  descriptors  in  Equation  2  are  poor,  as  seen  in  Table  4b.  Although  the 
correlation  and  cross-validated  correlation  coefficients  for  Equation  2  are  better  than  those  for 
Equation  1,  the  RMS  error  -  about  14°C  -  is  almost  three  times  larger.  Equations  1  and  2  are 
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also  similar  to  melting-point  QSARs  derived  for  other  ionic  liquids  in  that  they  involve 
hydrogen-bond  descriptors  and  NRI.9’17 

Ten  excellent,  single-descriptor  QSPRs  were  derived  for  the  densities  of  the  bromide  salts  and 
are  displayed  in  Table  5.  Most  descriptors  include  either  various  measures  of  hydrogen-bonding 
capability  or  have  to  do  with  the  partial  atomic  charge  distribution  and  are  interpreted  as 
characterizing  the  interactions  between  ions.  Another  descriptor,  the  average  electrophilic 
reactivity  index  for  a  carbon  atom,  is  similar  in  origin  to  the  previously  discussed  NRI19  and 
although  it,  too,  was  designed  to  characterize  a  chemical  reaction,  it  serves  here  to  quantify 
intereactions  between  ions.  Another  descriptor  -  a  principal  moment  of  inertia  normalized  by 
the  number  of  atoms  in  each  cation  -  is  a  measure  of  the  molecular  shape.  Each  descriptor 
predicts  very  accurately  the  liquid  densities  of  the  thirteen  bromide  salts  and  quantifies  very  well 
the  interactions  between  ions  while  in  the  liquid  state. 

Conclusions 

We  have  presented  QSPRs  for  the  melting  points  and  liquid  densities  of  a  new  class  of 
energetic  ionic  liquids  that  can  be  used  in  the  design  of  additional  l-substituted-4-amino- 1,2,4- 
triazolium  bromide  and  nitrate  salts..  Although  the  descriptors  used  were  originally  intended  to 
quantify  chemical  reactions,  we  believe  that  they  serve  here  to  describe  the  interactions  between 
ions.  Statistics,  while  certainly  more  than  adequate,  can  be  improved  by  designing  descriptors 
especially  for  ionic  liquids. 
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Table  1:  Melting  Points  and  Densities  of  Ionic  Liquids 


Substituent 

Bromide  Salt 

Nitrate  Salt 

Tm 

P 

Tm 

methyl 

92°C 

1.98  g/cmJ 

-54°C  (g) 

ethyl 

67 

1.69 

-55  (g) 

n -propyl 

63 

1.56 

33 

isopropyl 

92 

1.60 

66 

2-propenyl 

62 

1.59 

-50(g) 

n -butyl 

48 

1.46 

-50(g) 

methylcyclopropyl 

73 

1.58 

56 

//-pentyl 

54 

1.37 

29 

/? -hexyl 

76 

1.34 

0 

/7-heptyl 

94 

1.30 

35 

//-octyl 

80 

1.27 

34 

/7-nonyl 

81 

1.26 

53 

/7-decyl 

90 

1.23 

51 

Data  denoted  with  (g)  correspond  to  glass  transition  temperatures,  as  explained  in  the  text. 
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Table  2:  Comparison  of  Selected  Bond  Lengths  Found  in  Crystal  Structures  to  Calculated  Results 


H\  /H 
N4 


H 


N3 

/  \ 

C2  Cl' 

\  / 

N1 — N2 

/ 

C3 


-H 


Bromide  Salts 


ethyl  propyl  isopropyl  hexyl  heptyl 


Bond 

Crystal 

Calc. 

Crystal 

Calc. 

Crystal 

Calc. 

Crystal 

Calc. 

Crystal 

Calc. 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

N1-N2 

1.369 

1.349 

1.373 

1.349 

1.361 

1.349 

1.373 

1.349 

1.377 

1.349 

N1-C2 

1.312 

1.297 

1.307 

1.297 

1.298 

1.297 

1.308 

1.297 

1.314 

1.297 

N1-C3 

1.480 

1.478 

1.470 

1.477 

1.477 

1.492 

1.456 

1.478 

1.464 

1.478 

N2-C1 

1.310 

1.278 

1.307 

1.278 

1.297 

1.278 

1.293 

1.278 

1.303 

1.278 

N3-N4 

1.419 

1.383 

1.414 

1.383 

1.400 

1.383 

1.394 

1.383 

1.401 

1.383 

N3-C1 

1.356 

1.366 

1.364 

1.366 

1.355 

1.366 

1.350 

1.366 

1.357 

1.366 

N3-C2 

1.337 

1.327 

1.329 

1.327 

1.327 

1.328 

1.334 

1.328 

1.333 

1.328 
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Table  2,  continued 


Nitrate  Salts 


isopropyl 

methylcyclopropyl 

Crystal 

Calc. 

Crystal 

Calc. 

A 

A 

A 

A 

1.363 

1.349 

1.366 

1.349 

1.316 

1.297 

1.310 

1.297 

1.484 

1.478 

1.466 

1.477 

1.305 

1.278 

1.303 

1.278 

1.451 

1.383 

1.401 

1.383 

1.354 

1.366 

1.354 

1.366 

1.338 

1.327 

1.332 

1.327 
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Table  3:  Descriptor  Values  for  Melting-Point  QSPRs 


Substituent 

NRInh2 

HACA2 

Elumo 

HDCAlz* 

NRImm,c 

FHDCA 

au 

A- 

au 

methyl 

1.73  x  10'6 

4.95 

-0.876 

0.0161 

6.17  x  10'4 

0.208 

ethyl 

6.44  x  10'6 

4.29 

-0.991 

0.0144 

1.04  x  10'3 

0.220 

n -propyl 

3.13  x  10° 

4.35 

-1.049 

0.0120 

9.19  x  10'4 

0.215 

isopropyl 

1.10  x  10'5 

4.17 

-1.144 

0.0117 

1.27  x  10'3 

0.224 

2-propenyl 

1.83  x  10'5 

4.17 

-0.999 

0.0135 

1.17  x  10'4 

0.203 

n -butyl 

7.11  x  10'5 

4.39 

-1.09 

0.0109 

3.07  x  10‘3 

0.215 

methy  ley  c  lopropy  1 

6.17  x  I0'6 

4.06 

-1.11 

0.0119 

9.51  x  10'5 

0.215 

//-pentyl 

5.70  x  10° 

4.48 

-1.11 

9.97  x  I0'3 

1.21  x  10'3 

0.214 

//-hexyl 

2.61  x  10'5 

4.21 

-1.12 

9.37  x  I0'3 

3.11  x  10'4 

0.214 

//-heptyl 

1.28  x  10'5 

4.46 

-1.13 

8.60  x  I0'3 

1.16  x  10'4 

0.213 

/7-octyl 

6.72  x  10'6 

4.15 

-1.13 

8.50  x  I0'3 

5.30  x  10‘3 

0.217 

//-nonyl 

3.75  x  10'6 

4.10 

-1.14 

7.97  x  I0'3 

2.66  x  10'5 

0.217 

//-decyl 

2.19  x  10'6 

4.16 

-1.14 

7.61  x  I0'3 

1.46  x  10'5 

0.216 

*Zefirov's  empirical  atomic  partial  charges30’31  are  employed. 
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Table  4a:  Correlation  Matrix  for  Descriptors  in  Equation  1 : 


NRInh2 

HACA2 

1/Elumo 

NRInh2 

1.00 

- 

- 

HACA2 

0.168 

1.00 

- 

1/Elumo 

-0.110 

0.689 

1.00 

Table  4b:  Correlation  Matrix  for  Descriptors  in  Equation  2 


HDCAIha 

NRImin,C 

FHDCA 

HDCAIha 

1.00 

- 

- 

NRImin.r 

-0.359 

1.00 

- 

FHDCA 

-0.466 

0.205 

1.00 
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Table  5:  Density  QSPRs  for  l-substituted-4-amino-l,2,4-triazolium  bromides 


Intercept 

Coefficient 

r2 

F 

s2 

A 

Descriptor 

0.669 

3.24 

0.973 

397. 

0.00140 

0.954 

FHDSAZ:  Fractional  hydrogen-bond  donor  surface  area,  i.e.  hydrogen-bond 
donor  surface  area  divided  by  the  total  molecular  surface  area. 

0.974 

33.3 

0.972 

383. 

0.00140 

0.964 

ERIavg,c:  Average  electrophilic  reactivity  index  for  carbon  atoms.19 

0.756 

2.819 

0.964 

295. 

0.00180 

0.929 

FF1DSA:  Fractional  hydrogen-bond  donor  surface  area,  i.e.  hydrogen-bond  donor 
surface  area  divided  by  the  total  molecular  surface  area 

FF1ACA2:  sum  of  the  square  root  of  the  hydrogen-bond  acceptors’  surface  areas 

0.654 

75.4 

0.961 

270. 

0.00200 

0.951 

weighted  by  the  acceptors’  partial  charge.  The  sum  is  then  divided  by  the  total 
molecular  surface  area. 

16.3 

-16.0 

0.960 

262. 

0.00200 

0.922 

VH:  Average  valence  of  FI  atoms 

27.0 

-6.64 

0.945 

188. 

0.00280 

0.928 

V<>  Average  valence  of  C  atoms 

1.15 

48.4 

0.936 

162. 

0.00320 

0.915 

Ia/N:  Principal  moment  of  inertia  about  the  first  main  axis  divided  by  the 
number  of  atoms  in  the  molecule. 

2.06 

-1.26 

0.931 

148. 

0.00350 

0.906 

|ihC:Total  hybridization  component  of  the  molecular  dipole,  i.e.,  the  dipole 
calculated  from  partial  atomic  charges  subtracted  from  the  dipole  moment 
calculated  quantum  mechanically  in  Gaussian’03  as  the  first  derivative  of  the 
energy  with  respect  to  an  electric  field. 

FF1ASA2Z:  Fractional  hydrogen-bond  acceptor  surface  area,  i.e.,  the  sum  of  the 

0.744 

29.3 

0.924 

134. 

0.00380 

0.901 

hydrogen-bond  acceptor  surface  areas  divided  by  the  total  molecular  surface 

area.* 

0.727 

-87.2 

0.923 

132. 

0.00390 

0.892 

FNSA3Z:  Sum  of  the  negatively  charged  molecular  surface  area  weighted  by  the 
atomic  partial  charge.  The  total  molecular  surface  area  then  divides  the  sum. 

sR  "  "  "  "  OQ  O  1  ■  ■  ■ 

Zefirov's  empirical  atomic  partial  charges  ’  are  employed  as  weights  (FHASA2z  and  FNSA3z)  or  to  detennine  which  hydrogens 
are  donors  (FHDSAZ). 
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Table  6:  Descriptor  Values  for  Density  QSPRs 


Substituent 

FHDSAz 

ERIavgjc 

FHDSA 

FHACA2 

Vh 

Vc 

Ia/N 

Phc 

FHASA2Z 

FNSA3Z 

- 

au 

- 

A'1 

- 

- 

au 

au 

- 

- 

methyl 

0.390 

0.0298 

0.400 

0.0174 

0.898 

3.78 

0.0179 

0.0819 

0.0433 

-0.0138 

ethyl 

0.334 

0.0228 

0.356 

0.0134 

0.909 

3.81 

0.0116 

0.346 

0.0304 

-0.0105 

n -propyl 

0.269 

0.0183 

0.301 

0.0123 

0.917 

3.83 

9.14  x  10'3 

0.479 

0.0279 

-9.91  x  10'3 

isopropyl 

0.279 

0.0181 

0.311 

0.0120 

0.915 

3.83 

6.78  x  10'' 

0.426 

0.0258 

-8.77  x  10'3 

2-propenyl 

0.305 

0.0184 

0.305 

0.0125 

0.910 

3.85 

8.00  x  10'3 

0.347 

0.0287 

-0.0112 

n -butyl 

0.245 

0.0152 

0.239 

0.0115 

0.924 

3.85 

6.43  x  10'3 

0.500 

0.0265 

-8.68  x  10'3 

methyl- 

cyclopropyl 

0.267 

0.0152 

0.284 

0.0113 

0.915 

3.85 

6.84  x  10'3 

0.285 

0.0247 

-9.20  x  10'3 

//-pentyl 

0.225 

0.0130 

0.219 

0.0106 

0.928 

3.86 

5.48  x  10'3 

0.545 

0.0246 

-8.32  x  10'3 

//-hexyl 

0.211 

0.0114 

0.207 

9. 1 2  x  1 0'5 

0.932 

3.86 

4.20  x  10'3 

0.557 

0.0208 

-7.41  x  10'3 

n-heptyl 

0.195 

0.0101 

0.188 

8.87  x  10'3 

0.934 

3.87 

3.71  x  10'3 

0.594 

0.0208 

-6.85  x  10'3 

/7-octyl 

0.189 

9.10  x  I0'3 

0.185 

8.23  x  I0'3 

0.936 

3.88 

3.01  x  10'3 

0.608 

0.0186 

-6.19  x  10'3 

/?-nonyl 

0.177 

8.27  x  I0'3 

0.174 

7.62  x  I0'3 

0.938 

3.88 

2.72  x  10'3 

0.640 

0.0173 

-5.70  x  10'3 

/7-decyl 

0.170 

7.58  x  I0'3 

0.166 

7.25  x  10'3 

0.940 

3.89 

2.28  x  10'3 

0.655 

0.0166 

-5.57  x  10'3 
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Figure  1:  Molecular  structure  of  ionic  liquids  based  on  l-substituted-4-amino-l,2,4-triazolium. 
X-  is  either  Br~  or  NOr  and  R  is  a  series  of  aliphatics  or  other  common  organic  group  (seeTable 
1). 
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